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Recent experiments performed on inclusive electron scattering from nuclear targets have mea¬ 
sured the nucleon electromagnetic structure functions Fi{x,Q^), F 2 {x,Q^) and Fl{x,Q^) in 
^’^Al, and nuclei. The measurements have been done in the energy region of 

1 GeF^ < IF^ < 4 GeV^ and region of 0.5 GeF^ < < 4.5 GeF^ We have calculated 

nuclear medium effects in these structure functions arising due to the Fermi motion, binding energy, 
nucleon correlations, mesonic contributions from pion and rho mesons and shadowing effects. The 
calculations are performed in a local density approximation using a relativistic nucleon spectral 
function which includes nucleon correlations. The numerical results are compared with the recent 
experimental data from JLab and also with some earlier experiments. 

PACS numbers: 13.40.-f,21.65.-f,24.85.-|-p, 25.40.-h 

I. INTRODUCTION 

Charged lepton induced processes in the deep inelastic scattering region are used to probe quark and gluon structures 
of nucleons and nuclei. The inclusive lepton scattering from nucleon and nuclear targets is an important tool to study 
the nucleon structure and its modification in nuclear medium. The early results from the experiments performed at 
STAG in the kinematic region of high energy transfer(i/) and high momentum transfer((5^) corresponding to Deep 
Inelastic Scattering (DIS) exhibited remarkable phenomenon of Bjorken scaling [I|. In Bjorken scaling, the nucleon 
structure functions in the asymptotic region of very high are found to be independent of and depend only upon 
single dimensionless variable x{= 2 mt) instead of otherwise independent variables and v. This independence 
of nucleon structure functions led to the first evidence that nucleons consist of structureless constituents identified 
as quarks and gluons. Furthermore, the lepton-nucleon cross sections were found to be incoherent sum of elastic 
lepton scattering cross section from these structureless constituents confirming the asymptotic freedom predicted by 
QCD. At lower values of and x, the nucleon structure functions exhibit dependence which is attributed to the 
violation of Bjorken scaling due to higher twist correction in QCD M- These higher twist corrections arise due to 
quark-quark and quark-gluon interactions in the nucleon which play an important role when we move to lower 
from the asymptotic region in Q^. The nucleon structure functions in the scaling region of high Q^, and their 
evolution to lower region are described rather well using the methods of perturbative QCD and evolution equation 
of DGLAP Hi. 

In the case of deep inelastic lepton scattering from nuclear targets, the nuclear structure functions per nucleon are 
found to be quite different from the nucleon structure functions as discovered by the EMC effect first observed at CERN 
0 and confirmed by many experiments thereafter 0-0. These modifications in the nucleon structure function are 
due to the nuclear medium effects like Fermi motion, binding energy, nucleon correlations, mesonic contributions, etc. 
These are in addition to higher twist effects in QCD for nucleons and mesons in nuclear medium [sTI - fs^ . 

The nuclear medium effects in the DIS region are divided into four parts which are broadly identified in terms 
of Bjorken scaling variable x. These are (i) the shadowing effect which is effective in the low values of a;(< 0.1). 
It has been found that in this region R{=^) gets suppressed, and the suppression increases with the increase in 
the nucleon number A, (ii) the anti-shadowing effect which is effective in the region 0.1 < x < 0.3, where there is 
slight enhancement in the ratio i? which has been found to be independent of the nucleon number A. Shadowing and 
anti-shadowing effects are attributed respectively due to the constructive and destructive interference of amplitudes 
arising from the multiple scattering of quarks inside the nucleus, (hi) the EMC effect which is a large suppression in 
a wide range of a;(0.3 < x < 0.8), and is broadly understood as due to the modification of nucleon structure functions 
in nuclei, and (iv) the binding energy and the Fermi motion effect which is effective for x > 0.8 and this arises due 
to the fact that the nucleons in a nucleus are moving with an average momentum p < pp, where pf is the Fermi 
momentum. 


Corresponding author: zaidi.physics@gmail.com 



2 


In the DIS region, phenomenologically @ as well as theoretically various attempts have been made to 

understand the nuclear medium effects. Phenomenologically the studies have been made to obtain a nuclear correction 
factor by doing the analysis of the experimental data on charged lepton-nucleus scattering, (anti)neutrino nucleus 
scattering, pion-nucleus scattering, proton-nucleus scattering, Drell-Yan processes, etc. Theoretically many models 
have been proposed to study these effects on the basis of nuclear binding, nuclear medium modification including short 


range correlations in nuclei [39|-[63|, pion excess in nuclei |4lL 14,11.148l. I64l - l66| . multi-quark clusters [67H69| . dynamical 
rescaling lzO,[nl, nuclear shadowing [tI, [t^, etc. In spite of these efforts a comprehensive understanding of the 
nuclear modifications of nucleon structure functions valid for the entire region of x is still lacking . 

Recently JLab [l8j | has performed experiments using continuous electron beam facility with energies in the range of 
approximately 2-6 GeV, and precise measurements have been performed for the nucleon structure functions (x, Q^), 
F^{x,Q'^) and longitudinal structure function F^{x,Q'^) in the energy region of 1 GeV^ < W'^ < 4 GeV^ corre¬ 
sponding to low and moderate in the region of 0.5 GeV^ < < 4.5 GeV^ on several nuclear targets like 

^^Fe, etc. The modification of structure function in nuclear medium has also been studied earlier 

by SLAG [Tl|, NMC [l^, BCDMS [t^, etc. collaborations, in some of these nuclear targets as well as in a few other 
nuclear targets like etc. JLab also plans to upgrade electron beam energy to 12 GeV [zSIii and 

measure nucleon structure functions at low and moderate in the region of W relevant for the study of quark-hadron 
duality in nuclei. It is, therefore, important that nuclear medium effects are theoretically well understood in the deep 
inelastic region as well as in the resonance production region. 

Motivated by the recent experimental results [l8l| . we have performed a theoretical study to understand nuclear 
medium effect in nuclear structure functions F-f^j^{x, Q^), Fi^[x, Q^) as well as longitudinal structure function 
Fl ill nuclei for several nuclear targets like ^^Fe, ^^Gu, and in the region 

of moderate Q^. Our numerical results are compared with the available JLab data [l^ and also with the available 
results from earlier experiments performed by NMC [l^ collaboration. This study is a continuation of our previous 
study of nuclear medium effects in the extraction of electromagnetic and weak structure functions F^^{x,Q‘^) [^ . 
F^.^°'^{x, (3^)[5l,[5l,[63|, Paschos-Wolfenstein relation 61| and Parity Violating asymmetry in deep inelastic polarized 
electron scattering [6^ from nuclear targets. In the present paper, we have obtained the expressions of F^j^{x^Q'^) 
and F^ j~^{x,Q'^) separately without using the Callan-Gross relation. Therefore, these results can also provide a 
theoretical framework for testing the validity of Callan-Gross relation in nuclei. 

The present study has been performed with a microscopic model which uses relativistic nucleon spectral function 
to describe target nucleon momentum distribution incorporating Fermi motion, binding energy effects and nucleon 
correlations in a field theoretical model. The spectral function that describes the energy and momentum distribution 
of the nucleons in nuclei is obtained by using the Lehmann’s representation for the relativistic nucleon propagator 


and nuclear many body theory is used to calculate it for an interacting Fermi sea in nuclear matter 77[. A local 


density approximation is then applied to translate these results to a finite nucleus. Furthermore, we have considered 
the contributions of the pion and rho meson clouds in a many body field theoretical approach based on Refs. 

Due to the fact that JLab data [l^, [t^ have been taken in a region of relatively low we have not assumed the 
Bjorken limit. We have incorporated Target Mass Gorrection (TMG) following Ref. [^. This is effective at moderate 
and high x. The calculations are performed at the leading order(LO) as well as next to the leading order (NLO). 
The nucleon Parton Distribution Functions(PDFs) have been taken from the works of GTEQ group The NLO 
evolution of the deep inelastic structure functions has been taken from the works of Vermaseren et al. |8l| and van 
Neerven and Vogt (82l| . In the case of pions we have taken the pionic parton distribution functions given by Gluck 
et al. [1^ [S^. For the rho mesons, we have applied the same PDFs as for the pions as in Ref. EHl- We have also 
considered the effect of shadowing following Ref. 


M- 


The plan of presentation is the following. In Sect. [TTl some basic formalism for the inclusive lepton-nucleon scattering 
has been introduced. In Sect. IIIII the various nuclear medium effects have been discussed in brief, and in Sect. llVI the 
numerical results are presented. In Sect. El we conclude our findings. 


II. DEEP INELASTIC LEPTON-NUCLEON SCATTERING 


The double differential cross section for the reaction (depicted in FigH]) 

l{k) + N{p) ^ l{k') + X{p'), 1= e~, fj,~, 


in the Lab frame is written as 


dFlidE[ 


M T 

|k| 




( 1 ) 


( 2 ) 
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FIG. 1: Feynman diagram for the virtual photon induced deep inelastic lepton-nucleon scattering. l{l') is a lepton of four 
momenta fe(fe'), is a nucleon with four momentum p and X is the jet of hadrons with four momentum p'. 


where q = k — k' is the four momentum transfer, q^{= —AEiE'isin^^) < 0 , and a = e^/ 47 r. is the spin averaged 
leptonic tensor given by 




(3) 


where ^uE[ refer to the outgoing lepton. The hadronic tensor W'^'^ is defined in terms of , W 2 structure 
functions of the nucleon, 


= (14- - 


<7 




W, 


N 




P-9 I 
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M 2 


(4) 


Q 


with M as the mass of nucleon. 

In terms of the Bjorken variable x = 2^77 “ 2 M{Ei-e' 
transfer(= Ei — E[)^ the differential cross section is given by 


dxdy 


SMEiTra^ 


37 ^ and y = where = —q‘^ and v is the energy 
xyM \ 


xy^Fi\x,Q^)+[l-y- 


2Ei 


Ei\x,Q^) 


(5) 


where {x,Q^) = MWi{v^Q‘^) and F,^{x,Q'^) = are the dimensionless structure functions. In the 

Bjorken limit, i.e. —>■ 00 , e —>■ 00 , x finite, the structure functions Fjli 2(^1 Q^) depend only on the variable x and 
satisfy the Callan-Gross relation [s^ given by 2xF(^{x, (x, Q^). 

The total cross section for polarized photon (helicity A) interacting with unpolarized proton is expressed as 


crx = 


K 




( 6 ) 


where is the polarization vector of the photon with A = ±1,0; FT = ^ 2 A^ ^ ^ invariant energy of the 

virtual photon-proton system i.e VW^ = y^M'^ + 2M{Ei — E[) — Q"^ and W^i, is given by Eq.Q. Using the above 
expression, one may obtain the transverse and longitudinal cross sections as 

W^{„. f) 


(7t — 

(TL = 


K 

A'lr'^a 

K 


1 ± 




q^) - q^) 


(7) 


Now using the above expressions (i.e. Eq. Wi (v, q^) and W 2 {v, q^) may be written in terms of ctl and ctt, 
and the differential scattering cross section is written as 




dQ.idE[ 


Ft [crrix, ± eap^x, Q^)] = rTCTT(a;, [l ± eR'{x, Q^)] , 


( 8 ) 
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FIG. 2: (a) Lepton self energy (b) Photon self energy. The imaginary part is calculated by cutting along the horizontal line 
and applying the Cutkosky rules while putting the particle on mass shell. 


where F-r is the transverse virtual photon flux (= aE[/[A'K'^Q'^MEi{l — e)]), R' = e is the virtual photon 

polarization parameter given by 


e 


1 

l + 2(l + ^)tan2^ 


1 ,, 

t-2/-^ 


l-y + \ 




(9) 


^Lab is the Lab scattering angle. 

Using the above expressions for Ft and e, one may re-write Eq.Q as 

( 10 ) 

The expression for the longitudinal structure function E^(x, Q^) is defined as 

A l\/f 2^2 

Ff (:r, Q2) = (1 + Q^) - 2xE^{x, Q^) (11) 

where E(^{x, Q^) is purely transverse and E^(x, Q^) is a mixture of both. 

The nucleon structure functions are determined in terms of parton distribution functions for quarks and anti-quarks. 
In this work, for the nucleons we work at the leading order(LO) and next to leading order(NLO) and used CTEQ6.6 [s^ 
nucleon Parton Distribution Functions (PDFs). For evolution we have used the model of Refs. [nisi,113. It may 
be pointed out that we have also obtained the results at next to next to leading order i.e. NNLO (the results not 
shown here) and found the change between NLO and NNLO results to be less than half a percent. 




.N 


dElidE', 


= Fr 


dTT^a 


a;(lU2 _ m2) 


2xE^{x,Q^) 


(1 


2J2 


AM^x 


■)F, 


i^ix,Q^) - 2xF^{x,Q'^) 


III. NUCLEAR EFFECTS 


Now we evaluate the cross section for the reaction given by Eq. o for the lepton scattering taking place with a 
nucleon moving inside the nucleus. The expression of the cross section given in Eq.(I5]) is modified as 


d^idE'i 


o? |k'| 

7 M 


yiv 


1U7, 


( 12 ) 


where is the nuclear hadronic tensor defined in terms of nuclear hadronic structure functions W/^(i=l,2) as 


lU^'^ = 



FFi^ + 



PA-q 




PA-q A 


Ml 


(13) 


with Ma is the mass of nucleus. In the present formalism the lepton-nucleus DIS cross section is obtained in terms 
of lepton self energy in the nuclear medium. We evaluate the lepton self energy corresponding to the diagram shown 
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in Figl^l and the cross section for an element of volume dV in the nucleus is related to the probability per unit time 
(r) of the lepton interacting with the nucleons. TdtdS provides probability times a differential of area(dS) which is 
nothing but the cross section(dcr) i.e. 


da = TdtdS = r^dSdl = T-dV = r^dV = T^d^r, 
at V I k I I k I 


(14) 


where dl is the length of the interaction, v{= is the velocity of the incoming lepton and we have used k = vEi. 

Also probability per unit time of the lepton interacting with the nucleons in the medium to give the final state is 
related to the imaginary part of the lepton self energy i.e. 


Equating Eq. [H] and EqUSl one gets 


r 

2 


EiO^) 


ImS 


(15) 


Thus to get da for (/,/') scattering on the nucleus given by Eq. [T^ we are required to evaluate imaginary part of 
lepton self energy IrnS^k). For this, we write the lepton self energy (E(fc)) corresponding to the diagram shown in 
Fig 12^ by using Feynman rules 


— TT{k) = 


d Q N. n ■ ft -\- ni , / -s, x-rpo-/ N 

-j u/(k) iej>^ *67 ui{k) —^ (-z) {q) - 


(27r) 


q 


(17) 


which for unpolarized leptons can be written as 


E(fc) = ie^ j 


d‘^q 


1 1 

(27r)^ q'^ 2m 


1 


/fc'2 _ 


+ ie 


n^"(g), 


(18) 


where T\^'^{q) the photon self energy and is the leptonic tensor given by Eq|3l Now we shall use the imaginary 
part of the lepton self energy i.e. ImT,{k), to obtain the results for the cross section and for this we apply Cutkosky 
rules liMl 


E(fc) —>■ 2i ImT,{k) 

D{k') —>• 2i6{k'°) ImD{k') (boson propagator) 

\T^^’'{q) 2ie{q°) ImUf^-^iq) 

G{p) —>■ 2i9{p°) ImG{p) (fermion propagator) 

which leads to 


( 19 ) 


( 20 ) 


Notice from Eq. [501 E(fc) contains photon self energy 11^'^, which is written in terms of nucleon propagator Gi and 


meson propagator Dj and using Feynman rules [88| following Fig. [2] this is given by 


n^'^(g) 


d'^p 


G(p)E En,_, 


d^p'^ 

(2^)4 


i[Gi{p[) n D,{p'^) 

I 3 


N 

< X\J^\H X X| J"|i7 >* (27r)4 5\q+p-Y,p[), 


( 21 ) 


where Sp is the spin of the nucleon, Si is the spin of the fermions in X, < X| J^|iJ > is the hadronic current for 
the initial state nucleon to the final state hadrons, index I runs for fermions and index j runs for bosons in the final 
hadron state X. 
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FIG. 3: Nucleon self-energy in the nuclear medium 


To derive the nucleon propagator in nuclear medium G(p) we start with the relativistic free nucleon Dirac propagator 
G^{Pq,p) and write it in terms of the contribution from the positive and negative energy components of the nucleon 
described by the Dirac spinors u(p) and u(p) using their appropriate normalization (uu = 1), which results [s^: 


M j' Sr^r(p)Mr(p) ^ Sr 'Cr (-p)fr-("P) j 
E{p) \po - E{p)+ie po + E{p)-ie J 


( 22 ) 


where M is the mass and i?(p) is the relativistic energy(A/p2 -|- M^) of an on shell nucleon. We shall retain only the 
positive energy contributions as the negative energy contributions are suppressed. 

Moreover, the relativistic nucleon propagator in a non-interacting Fermi sea may be written as 


G°(po,p) 


M 

W) 



1 - »(p) 
po - E{p) + ie 


njp) 

Po - E{p) - ie 


+ 


Y.r'^r{-P)vr{-P) 
Po + E{p) - ie 


(23) 


where u(p) is the occupation number of the nucleons in the Fermi sea, n(p)=l for p< pp^ while n(p) =0 for p> 

In the interacting Fermi sea, the relativistic nucleon propagator is written using Dyson series expansion (depicted 
in Fig Bt in terms of nucleon self energy E^(po, p)- This perturbative expansion is summed in a ladder approximation 
to giveQ: 


G{p) 


M 'Er'^r(p)Ur(p) M Y.r Urip)Urip) ^ ^ Es ««(P)'Ws (P) 

E(p) po-E{p) Eip) po-E(p) ^^°’P^i?(p) po-E(p) 

M Ur(p)Mr(p) 

^^(P) ^ PO- E(p) - Ur(p)^^ipo, P)Mr(p) 


(24) 


One may notice from the expression for the nucleon propagator G(p) given in Eg 1241 that it contains nucleon self 
energy Yl^(po,p). The nucleon self-energy is written using the techniques of the standard Many-Body Theory [za. 
The inputs required for the NN interaction are incorporated by relating them to the experimental elastic NN cross 
section. Furthermore, RPA-correlation effect is taken into account using the spin-isospin effective interaction as the 
dominating part of the particle-hole (ph) interaction. Using the modified expression for the nucleon self energy, the 
imaginary part of it is obtained. Due to RPA effect the imaginary part of the nucleon self-energy is quenched specially 
at low energies and high densities and this also depends on nucleon energy po as well as nucleon momentum p in the 
interacting Fermi sea. The imaginary part of the nucleon self-energy fulfills the low-density theorem. The real part of 
the nucleon self energy is obtained by means of dispersion relations using the expressions of the imaginary part. The 
Hartree and Fock pieces of the self-energy do not contribute to the imaginary part and therefore cannot be obtained 
by means of the dispersion relations and are explicitly added. For the real part, the only added piece is the Fock term 
assuming that the interaction contains spin-isospin excitations having longitudinal and transverse components. This 
model, of course, misses some pieces of Hartree type, which depend on the density but not on po or p, but its effect 
is expected to be small in the study of present interest. Following Ref. [i^, the relativistic nucleon propagator G(p) 
in a nuclear medium is thus expressed as: 


G(p) = 


M 

E{p) 


^Mr(p)Ur(p) 



du! 


Sh{uJ,p) 
Po-uj -ir] 



Sp{u},p) 
acu - 

Po - cj 


(25) 
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FIG. 4: Sh{ijJ, p) vs u for p < pj?(Left panel) and p > pi?(Right panel) in ^^(7(801^ line) and ®®f’e(daslied line). 


where Sh{u!,p) and 5^,(0;, p) being the hole and particle spectral functions respectively, which are given by [48 


Sh{po,p) = 


Mv) 


ImS^(po,p) 


{po-E{p) - ^ReS^(po,p)) + (^ImS^(po,p)) 


(26) 


for po < p. 


Sp{po,p) 


TT 


4|yImS^(Po,p) 

(po-E{p) - ^ReS^(po,p)) + (^ImE^(po,p)) 


(27) 


for pq > p. p is the chemical potential given hy p = + ReT,^ present numerical calculations, 

the expressions for which have been taken from Ref. ffl\ . A few properties related with spectral function have been 
presented in Appendix-El In FigHl following Ref. [i^, we have shown Sh{oj, p) vs co (where w = po — M), for p < pp 
and p > pf in and nuclei. It may be observed that for p < pp the hole spectral function Sh almost mimics 
a delta function as it corresponds to a Lorentzian distribution with a very narrow width. While for p > pp, Sh is not 
exactly zero, although very small in magnitude but has a longer range. This behaviour is different from independent 
particle model where it is exactly zero and this difference arises due to nucleon correlation [s^. 

Using Eas ll6l and 1201 and performing the momentum space integration one may write the cross section as: 


da^ 


dflidE', 




(28) 


At this stage if Ea. dT^ and Eg. (1^5)) are compared it may be inferred that is related with as 




ImIE‘''d^r 


(29) 


Using Ea. (l2^ and the expressions for the free nucleon and meson propagators in Ea. d^ . and finally substituting 
them in Ea. (l29L we obtain nuclear hadronic tensor for an isospin symmetric nucleus in terms of nucleonic hadronic 
tensor and spectral function, given by 


Wf = A 


d-^r 


d^p M 


(27r)3 £;(p) 


dpoSh{po,p,p{r))W^^{p,q), . 


(30) 

















In this way we have incorporated Fermi motion, Pauli blocking and nucleon correlations through the inclusion of 
spectral function Sh{poiP, pif))- The normalization of the spectral function and the quantities obtained from it are 
given in Appendix-lAl 

Accordingly the dimensionless nuclear structure functions 2 ( 2 ;, are defined in terms of 3(^7 0^) as 


F,^{x,Q^) = Ma 

where 


^ = 90 j P^a~ ^ad is the mass of a nucleus. 


Ma 

Taking the xx component of Eq.(j4]), we have 


IF. 


N 


_ f qxQx 


IF, 


N 


^ ( P-^ \ ( P-^ \ wiV 

M2 q2 J ^2 ) ^2 


Choosing q along the z-axis, we obtain for = (go, 0, 0, q) and = {E^, p), 

IFiI(x^, Q2) = wf{xN,Q^) + -^plwi^ixN, Q") 
Similarly taking the xx component of Eq. (US and using p^ = {M^ ,0), we have 

Wi{XA,Q^) = W^{xA,Q^) = 

Using equations Eg. (1551) and Eg. (1551) in Eg. (1501) . we have 


F,^{xa,Q^) 

AM 




M2 


Px 


(2^)3 E(p) 

F^{xn,Q^)' 


dpoSh{po,P, p{r)) 


rEf(crjv,g2) 


M 


where xjsi = 


_ Q _ 


Q" 


and xa = ^ - — 


2p-q 2(paqo-Pzqz) A A A 2AIqo ' 

Eor nonisoscalar nuclear target the above equation is written as 


dpoSlipo^p.p'^ir)) 


T—p,n 

, 1 „ 2TI(a;jv,g^] 

+ M2^" 


F[iXN,Q^ 

M 


Similarly taking the zz component of Eq.(|4l) and Eg. (115)) . we obtain 


iNn,2\ 9o wAf/™ ni2\ , ^ ((Pz^ PA ^z) \ ^2\ 


IFf,(crw,Q^) = ^W,^ixN,Q^) +j^ 


W2^ixN,Q^) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


IF/,(x^,Q^) = :^W^\(xa,Q^) +^W^^Jxa,Q'^) 


Using equations Eg. (1571) and Eg. (1551) in Eg. (1501) . we have 

W^\ixA,Q^) +^W^\{xa,Q^) = 2 ^ f d^r f 
d ^^ J J 


9o 


d?p M /■'' 


(27r)3 E(p) J_^ 
2 


dpoShipo,p,p"'{r)) X 


(38) 


T—p,n 

'"^o TT.r^ ^ 2 ^ , 1 ( {Pzq^ - P-q dz) \ ^ 2 ^ 

^ IFi {XN, Q ) +J^ 1 - -A - ) ^2 {XN, Q ) 


( 39 ) 
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fl'l'l 


FIG. 5: Lepton self energy diagram including particle-hole(lp-lli), delta-liole(lA-lli), Iplh-lAlh, etc. 


excitations. 


With the expressions of and in terms of F^j^{x,Q'^) and F^j^{x,Q'^), we obtain from 

Eg. (15^ following expression for F^j^{x, Q^): 


= 21 : 


(SF ^ /I 


\Q^ 

- pI\ 

.ql 

^ 2M2 ) 


(po - Pz lY 

M 2 


Pz 


{po - Pz i)qoqz 


+ 1 


M 


Po - Pzl 




(40) 


where 7 = 

' qo _ _ 

Thus for the numerical calculations we shall use Eg. (1551) for Ff'^{x,Q'^) and Eo. (|5(I)) for F^^{x,Q'^). These 
numerical results will be labeled as the results with spectral functions. Now we include the mesonic contributions 
coming from pion and rho mesons. We follow the same approach as in the nucleon case. However, for completeness, 
we explicitly present a brief description of this formalism. 


A. Mesonic contributions 

There are virtual mesons associated with each nucleon bound inside the nucleus. These meson clouds get strength¬ 
ened by the strong attractive nature of nucleon-nucleon interactions. This leads to an increase in the interaction 
probability of virtual photons with the meson cloud. The effect of meson cloud is more pronounced in heavier nuclear 
targets and dominate in the intermediate region of x(0.2<x<0.6) which leads to an enhancement of nuclear structure 
function F ^2 [13 ■ To obtain the contribution from the virtual mesons, we again evaluate lepton self energy and for 
this a diagram similar to the one shown in Fig|2]is drawn, except that instead of a nucleon now there is a meson which 
results in the change of a nucleon propagator by a meson propagator. This meson propagator does not correspond 
to the free mesons as one lepton (either electron or muon) can not decay into another lepton, one pion and X but 
corresponds to the mesons arising due to the nuclear medium effects by using a modified meson propagator. These 
mesons are arising in the nuclear medium through particle-hole(lp-lh), delta-hole(lA-lh), Iplh-lAlh, 2p-2h, etc. 
interactions as depicted in FigIS] 

In the present model, we have considered the contribution from tt and p meson clouds following the many body field 
theoretical approach as used in the case of bound nucleons. We shall make use of the imaginary part of the meson 
propagators instead of spectral function. In the case of pion following Ref. [13, we replace 

-27i--^^S'/i(po,p) W ^^{ p , q ) 
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FIG. 6: F^{x,Q'^) vs X at a fixed Q^, in A='^^C, ^^Fe, and with spectral function(dashed-dotted line) and full 

calculation (dashed line) at LO and the results with spectral function(dashed-double dotted line) and full calculation(solid line) 
at NLO. Experimental points are the JLab data [T^ . 


in Eg. (HOI) by 

ImD{p) 0{po) 2IE“^(p,q) 

where D{p) is the pion propagator in the nuclear medium given by 

Dip) = [po'^ - - ml - n^(po,p)]"^: 


(41) 


with 


P/mlF\p)p^n* 

1 - p/mlV[n* 


(42) 


Here, E(p) = (A^ — mJ)/(A^ + P^) is the ttNN form factor, A=1 GeV, f = 1.01, V[ is the longitudinal part of the 
spin-isospin interaction and H* is the irreducible pion self energy that contains the contribution of particle - hole and 
delta - hole excitations. 

Following a similar procedure, as done in the case of nucleon, the contribution of the pions to hadronic tensor in 
the nuclear medium may be written as 

WZ =^JdP J ^ eipo)i-2) ImDip) 2m^wp{p, q) (43) 

A factor of 3 arises due to the three charge state of pions and a factor of 2 is absent as compared to the nucleon as 
the pions are spinless particles. 

Ea. (H51) also contains the contribution of the pionic contents of the nucleon. Since these pionic contents are already 
contained in the sea contribution of nucleon, therefore, the pionic contribution of the nucleon is to be subtracted 
from Ea. (l43|) . in order to calculate the contribution from the pion excess in the nuclear medium. This is obtained by 
replacing ImDip) by 6ImDip) [i^ as 


ImDip) —>■ SlmDip) = ImDip) ~ P 


d ImDip) 
dp 


lp=0 


(44) 
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FIG. 7: Results of full calculation at LO(LO Total) and NLO(NLO Total) for 2a:F’]^(a:, Q^)(dashed-dotted line: LO Total and 
solid line: NLO Total) for A=^^C', are presented at different Q^. Experimental points are the JLab 

data [T^. 


which leads to 


F^^{xTr,Q‘^) = -QAM J (fr J SImD{p) 2m^ 


Fl-jz {X-K J Q ) 


TTItt 


IpP - pI 

2(po go - Pzqz) m-rr 


(45) 


where 


Following the same procedure as for F^j^{x, Q'^) obtained in Ea. (l40l) . the expression for F^^{x), is given by 

F^^{x^,Q'^) = -6 J d^r J -^^e{po) 61 mD{p) 2m^ x 



hpp - pn 

1- 

1 

bf to 

to 

^lO 


+ 


Pz 


{po - Pz i)qmz 


+ 1 


Po - Pz 1 


F2Ax.,Q^) (46) 
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j^Modified / ipSF/ /^2\ 

FIG. 8; Vi = —^- Modified, — 2 ' - (i=li2) in % vs x, at = 2 GeV^ in ^^Fe. Here Ff^ stands for nuclear structure 

functions F^{x,Q^) {i = 1,2) obtained using spectral function only and stands for nuclear structure functions 

Ff{x,Q^) {i = 1,2) evaluated (i) with mesonic effects along with the spectral function, and (ii) when shadowing is also 
included in (i). The solid dotted(solid) line is the result obtained for n using case i(ii) and dashed dotted(dashed) line is the 
result obtained for r 2 using case i(ii). 


Similarly the contribution of the p-meson cloud to the structure function is taken into account in analogy with the 
above model and the rho structure function is written as [d^ 


= -I22IM J d^r J SlmDpip) 2m^ 


Fip{xp,Q‘^) 


- pI F2p{xp,Q'^) 


rrio 


2(po 9o - Pz dz) rup 


(47) 


FtM^) 


/ /» 

d^r J 5ImDp{p) 2mp x 



(\p\^ - pI) 

, {Po - Pz if 

( PzQ^ 

d^z 

1 2m2 

+ 

V (po - Pz i)qoqz 



rrip 

Po - Pzl 


F2p{xp,Q'^) 


(48) 


where Xp = and Dp{p) is now the p-meson propagator in the medium given by: 

Fp{p) = bo^ - - TOp - n*(po,p)]"^, 

where 

P/mlCpF^^{p)p^n* 

^ 1 — /^/rnpH^n* 


(49) 


(50) 


Here, is the transverse part of the spin-isospin interaction, Cp = 3.94, Fp{p) = (A^ — m^)/{A^ + P^) is the pNN 
form factor, Ap=l GeV, / = 1.01, and H* is the irreducible rho self energy that contains the contribution of particle 
- hole and delta - hole excitations. Quark and antiquark PDFs for pions have been taken from the parameterization 
given by Gluck et al.jl^ and for the rho mesons we have taken the same PDFs as for the pions. 

We have tested that our model fulfills the momentum sum rule as expressed in Eq.(88) of Ref (including also 
the p meson). The procedure is straight forward by using Eqs. 36-37 from Kulagin and Petti and Eq.(lO) from our 
paper published in Ref. [^ . The pion y,, and nucleon fractions of the light cone momentum are related by 


Vt^+Vn 


Ma 

AM' 


(51) 
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FIG. 9: Results of the full calculation at LO(LO Total) and NLO(NLO Total) for the longitudinal structure function F^{x,Q^) 
vs X, in various mic\ei{A=^^C, Al, and ^^Cu). These results are presented at different Q^. Experimental points are 

JLab data [3. At NLO, these results are also presented with the spectral function(NLO SF). Nuclear targets ^^Fe and 
are treated as nonisoscalar in the present calculation. 


where Ma is the nucleus mass. The nucleon quantities can be easily obtained from the spectral function. Our results 
are like this: For iron yjv = 0.967, 7r+p should account for 0.024. 


IV. RESULTS AND DISCUSSION 

For the numerical calculations we have used Ea. dMl) for Q^) and Ea. (l40l) for Q^). The contributions 

of TT and p mesons have been taken into account using Eg. lHSl) for F^^{x,Q‘^), Ea. (|T5)) for F.^^{x,Q‘^), Eq. dTTl) for 
Fip{x,Q'^) and Eo. lHSl) for Q^). We have added the mesonic contributions to the nucleonic contribution 

to get the full result. It should be mentioned that all the parameters entering in the above mentioned equations 
have been fixed by our earlier works and have not been taken as free parameters. We have used CTEQ 

parameterization [80| for the quar k/antiquark PDFs for nucleon and the quark/antiquark parameterization for pion 
from the works of Gluck et al. [^. For the rho meson we have taken the same PDFs as for the pion. Moreover, while 
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determining nucleon structure functions in terms of parton distribution functions, it is important to include target 
mass correction(TMC) to take into account the effects associated with the non-zero mass of the target. TMC effects 
are more pronounced at large x and moderate Q^. We have taken into account the target mass correction following 
the works of Schienbein et al. . The shadowing effect has been incorporated following Ref. [5l| and is pronounced 
at low x(x<0.1). Therefore, in the present studied region of x(> 0.1) the effect of shadowing is almost negligible and 
the results for this are not explicitly shown in all the figures. 

The results are presented for different cases. The first case is when we take the contribution of nucleon spectral 
function only i.e. Eg. (1551) for and Ea. (li(l for which includes the effect of Fermi motion, 

nuclear binding and nucleon correlations, and perform the calculation at the leading order(LO). The numerical results 
are referred as LO(SF). We then include the contribution from meson clouds as well as shadowing effect and this we 
call as the results of our full model. This we refer as LO(Total). These results have also been obtained at NLO for 
both the cases. The results in these cases are referred as NLO(SF) and NLO(Total), respectively. 

In Figini we have presented the results for F^{x,Q^) in various nuclear targets like and 

at different Q^. We find that in at = 2 GeV^, the results obtained using full model at LO, is about 18% 
higher at x=0.2 in comparison to the results obtained using spectral function only. This difference decreases with the 
increase in x, for example at x=0.4, it is 12% and becomes almost negligible at x=0.6. We observe that this difference 
increases with the increase in mass number i.e. for it is ^ 20% higher at x=0.2 and 14% at x=0.4, while in 
it is ^ 23% higher at x=0.2 and 17% at x=0.4. In ^^®S'n(not shown in the figure) this difference becomes about 25% 
at x=0.2 and 18% at x=0.4, whereas in (not shown in the figure) this difference increases to 26% at x=0.2 and 

20% at x=0.4. The difference is mainly due to mesonic effects and is negligible for x > 0.6 for the nuclei considered 
here. When the results obtained by using the full model at NLO are compared with the results evaluated at LO, we 
find that the results decrease from the LO values. For example, in the case of it is lower by 10% at x=0.2 and 5% 
at x=0.4. Then there is a cross over around x=0.5 and beyond that it increases with increase in x. For example, it is 
around 35-38% larger at x=0.8 for the nuclei considered here. The effect of shadowing has been found to be almost 
negligible in the presently studied region of x. For example, the effect of shadowing is around 1-2% at x^O.2 in 
which increases to 2-3% for heavier nuclei like and We have also shown in this figure, JLab experimental 

data [l^ and find that our results obtained with full model at NLO agree reasonably well with the JLab data. 

In Figl3 we present the results for 2xF^{x,Q‘^) obtained using Eas. (l36t . (|45|) and (iTTt . The results are presented 
with full model at LO and NLO in several nuclei like ^"^Al, ^^Fe and at various values of Q^. We find 
that in general these results are qualitatively similar in nature to that found in the case of F^{x,Q^), however, 
quantitatively there is some variation, specially in the region of low x where mesonic effects play a role. The results 
are also compared with the JLab data 18| and we find a reasonable agreement with the experimental data. 

We separately show the effect of mesonic contribution and shadowing effect in F^{x,Q‘^) and F^{x,Q^), by pre- 

senting the results for n = - Modified,- —^i(i=l!2) in FiglU These results are shown in where 

Ff^{x,Q^) stands for the results obtained for the nuclear structure functions using the spectral function only while 
pModified is the result obtained when we include (i) mesonic(7r -|- p) contributions and (ii) mesonic(7r -|- p) 
contributions and shadowing effects. At low we find that the mesonic contributions are larger which become 
smaller with the increase in . For example, at = 2 GeV^, the mesonic contribution is 24-28% at x=0.2, which 
becomes 2-4% at x=0.6. Similarly at = 5 GeV^(not shown in the figure), the mesonic contribution is found to 
be 20-24% at x=0.2, which becomes 1-2% at x=0.6. We also find that the mesonic contributions to F^{x,Q^) at a 
given x and is larger than the mesonic contribution in F^{x, Q^) over the whole range of x and Q^. The results 
presented in Figsin]and[7]for F^{x, Q^) and 2xF^{x, Q^) give useful predictions for F^ 2 i^^ which may 

be tested in future experiments. It may be observed that the mesonic contributions are really different for F^{x, Q^) 
and F.^{x,Q^) structure functions and when shadowing contributions are included both curves get reduced because 
the shadowing goes in the opposite direction to the enhancement due to meson cloud contribution and this can be 
viewed as a kind of compensating effect. 

In Figini we present the results for the longitudinal structure function F^{x, Q^) using the full model at LO as well 
as at NLO. Here we have also shown the results with the spectral function only, obtained at NLO. These results are 
presented for various nuclear targets like and at different Q^. We find that mesonic contributions 

to F^{x, Q^) are quite large at low x and and becomes negligible for x > 0.6 for almost all values of studied 
here. For example, in ^^G at = 2 GeV^, the mesonic contributions is ^ 58% (of total F^{x, Q^)) at x=0.2 which 
reduces to 20% at x=0.5. The mesonic contributions increases for all x(x < 0.6) with the increase in mass number as 
we go from to ^^Gu as shown in FigI3 Thus the results presented in this figure are also an equivalent proof of 
the fact that the mesonic contributions to F^{x, Q^) and F^{x, Q^) are really different. This is because, in principle, 
F^{x,Q^) should be zero if Callan-Gross relationship were exactly fulfilled. While the nucleonic contributions(SF 
only) of F 2 and 2a:Fi show a more prominent trend to cancel while computing F^, whereas the mesonic contributions 
when added to it leads to a significant behaviorial change which results in a broken Callan-Gross relation(compare 
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FIG. 10: Results obtained using the spectral function(NLO SF) and the full model(NLO Total) at NLO for the ratio 
^wc}^’n 2 \ i 1=1,2. Experimental points are the NMC result [l^ . 


the curve shown by the dotted lines to the curve shown by the solid line in Fig[^. The results at NLO are larger than 
LO for X > 0.4 and are in better agreement with the experimental results for F^{x, Q^) for JLab [l^. However, there 
are some exceptions specially at low (< 2 GeV'^) as seen in the case of ^^Al and ®®Fe. At smaller in the region 
of < 2 GeV^, the non-perturbative QCD effects and their possible enhancement in nuclear medium may play an 
important role which are beyond the scope of this work. However the present work makes important prediction for 
F^{x,Q^) in lower x region which may be tested in future experiments. 

We present the results for and in FiglTOl using the spectral function and the full model at NLO 

and compared the results with NMC data 13 . We find that the results are in better agreement with the experimental 
observations when mesonic contributions are included in addition to the nuclear medium effects. We find a strong 
nuclear dependence as the absolute values of the slope increases with the increase in the mass number. This may be 
noted that although mesonic cloud contributions seem to be very different for F^{x,Q'^) and F^{x,Q^) separately, 
they also seem to cancel when computing the ratios between structure functions for different nuclei. 


V. CONCLUSION 


To conclude, we have studied nuclear medium effects in electromagnetic nuclear structure functions Ff-{x, Q^) and 
F^{x, Q^), and the longitudinal structure function Fp{x, Q^). For the nuclear medium effects, we took into account 
Fermi motion, nuclear binding, nucleon correlations, effect of meson degrees of freedom, etc. The calculations are 
performed both at LO and NLO. 

The theoretical expressions for F^{x, and F.f{x, Q^) have been obtained without assuming Callan-Gross relation 
at nuclear level but it has been assumed at nucleonic(and mesonic) level while computing Fi 2 (x^ Q^) and F^ 2 ^{x, Q^) 
for free nucleons and mesons. The theoretical results are then compared with the experimental data for F^{x,Q^), 
F^{x,Q'^), Fp{x,Q^) from JLab data [l8| and found in good agreement with them (for x > 0.4) except for Fp at 
very low < 2 GeV^. The present work makes predictions for nucleon electromagnetic structure functions at lower 
X in various nuclei in the region of 2 < < 4GeU^, which will be useful in analyzing the future experiments being 

done for studying the nuclear medium effects in lepton nucleus scattering at low and moderate Q^. The results are 

also compared with the old NMC data on JyA with good agreement and predictions are made for other nuclei. 
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Appendix A: Properties of the spectral function 


The hole and particle spectral functions fulfill the following relations, 

f dpo Sh{po,p) = n{p) 


(Al) 



dpo Sp{po,p) = 1 - n(p) 


(A2) 


where n(p) is the Fermi occupation number. 

Therefore, the following sum rule is also full filled 



dpo Sh{po,p) + / dpo Sp{po,p) = 1 


(A3) 


In the absence of interactions, the nucleon energy po is the free relativistic energy E{p) and the dressed propagator 
G{p) reduces to the free propagator G°(p) i.e. if = 0 then 


Sh{po,p) = Sp{po,p) = S{po - E{p)) 


then 


J dpo Sh{po,p) = J dpo <5(po - i^(p)) = I Q |^<^(pj 
J dpo Sp{po,p) = J dpo (5(po - A(p)) = I Q m>F(p) 


(A4) 


(AS) 

(A6) 


Thus in the limiting case of vanishing self energy the expression for the spectral functions given in Fasl^and[?7l 
collapse to a representation of Dirac delta function. If A(p) is the total relativistic energy, then p must have the 
nucleon mass M incorporated i.e. 

p = M + ep (A7) 

With this definition, we can perform a constant shift in the integration variable po, given by: 

Po = ui + M (A8) 


And with this shift, the integrals stand: 

p — M > E{p) — M ^ ep > e(p) 
p — M < eIp) — M ^ ep < e(p) 

p — M < E{p) — M ^ ep < e(p) 
p — M > eIp) — M ^ ep > e(p) 

where e(p) = Eijp) — M is the nucleon kinetic energy, which in the non-relativistic regime can be approximated by 


rP 

pfi—M 

f 1 if 

/ dpoSh{po,p) = 

' —OO 

/ duj 6{oj + M — E{p)) = 

J —OO 

{ 0 if 

^OO 

/•OO 1 

' 1 if 

/ dpo Sp{po,p) = 

J ^ 

/ dw S{uj + M — E{p)) = < 

Jp-M \ 

0 if 


e(p) « 


P 

2M 


(A9) 
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If ef = is the Fermi energy, then the two step functions 0{pF — \p\) and 6 *(|p| — pf) are the solutions of the above 


2M 

integrals, namely 


f dpo Sh{po,p) = e{pF - \p\) = no{p) 

J —oo 

pOO 

/ dpo Sp{po,p) = 6 »(|p| -pf) = 1 - no(p) 
J U 


(AlO) 

(All) 


Thus in the absence of interactions, the full dressed propagator reduces to free one and one may write 

M 


Gipo.p) = 


Eip) 

M 

W) 

M 


^Mr(p)u^(p) 

r 

^M^(p)u,.(p) 




' —OO 


- U - IT] 

5{UJ- £;(p)) 


Ip P0-UJ + iri 


/ du) 

1-00 Po-w-ir] 


^Mr(p)Ur(p) 
r 

£;(p) 


Eip) 

M 


dw 


IpI -Pf) 


5{uj- £;(p)) 

Po - w + Z77 


^ipF - IpI) _ 

Po - E{p) -ir] po - E{p) + ir] 

»o(p) 1 - 'rto(p) 

Po - E{p) - ir] Po- E{p) + ir]_ 


(A12) 


The hole spectral function Sh{po,p) is physically interpreted as equal to the joint probability of (i) removing of a 
nucleon with momentum p from the correlated ground state, and (ii) of finding the resulting system of (A- 1 ) nucleons 
with an energy in the interval po and po + dpo. 

The normalization of this spectral function is obtained by imposing the baryon number conservation following 
Ref. ili: 


and 


(7V|R'"|Af) = m(p) 7 ^u(p) = R = 1 , p^ = (£'(p),p) 


{A\B^\A) = - I ^VrTr[G{po,p)r]e 


(A13) 


(A14) 


where V is the volume of the normalization box and expliporj), with 77 —>■ 0'*', is the convergence factor for loops 
appearing at equal times. 

Using the expression given by Ea. ipS]) . it can be seen that the convergence factor limits the contribution to the hole 
spectral function and one gets 


{A\Bf^\A) = y j 


= 2V 


d?p M 
(2f)3 E{p) 

f d^p M 


Tr 


^Ur{p)Ur{p)-/^ 


-Tr 




(27r)3 E{p) 

(27r)3 E{p) M 


2M 

■A* 


-T 


f Shiuj,p)duj = B 

J —OO 


Sh{uj,p)duj 


Shiuj,p)du; 


- u Pa 


Ma 


(A15) 


It is to be noted that in the last step we have imposed that this matrix element gives the right current with B baryons, 
in analogy to the expression given by Ea. (IA13l) . p(^ is the momentum of the nucleus. The operator {p + M) comes 
from ^ Mj.(p)'Ur(p) which depends only on p, and that corresponds to a free particle with p'^ = (E(p), p), therefore, 
the operator (^ + M) when the particle is on shell is written as 




(A16) 


We have ensured that the spectral function is properly normalized and checked it by obtaining the correct baryon 
number and binding energy for a given nucleus. In the local density approximation, the spectral functions of protons 
and neutrons are the function of local Fermi momentum. The equivalent normalization to Eq. (IA16I) is written as 


f d?p 

J JEY j = PpA'^) 


(A17) 
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Nucleus 

<T> jA (MeV) 

B.E./A (MeV) 


20.2 

7.6 

"Ml 

27.8 

8.3 

^'^Fe 

30.0 

00 

bo 

'^^Cu 

29.3 

8.7 

"“S'n 

31.8 

8.6 


33.7 

7.9 


32.7 

7.8 


TABLE I: Kinetic energy per nudeon(< T > /A) and binding energy per nucleon/A) for Al, ^^Fe, ^^Cu, 

and 


where factor 2 is due to two possible projections of spin i particle. is the Fermi momentum of proton(neutron) 

1 / Q 

inside the nucleus which is expressed in terms of proton(neutron) densities given by = [37r^/9p(„)(r)] 

These nucleon densities are in turn related with the nuclear densities /9(r)(like Pp{r) = ^p{r) and p„(r) = 
the parameters of which are determined from lepton scattering experiments. In the present calculation we have used 
harmonic oscillator density for nucleus and two parameter Fermi density for and 

nuclei which are taken from Refs. [9l|-[9^. This leads to the normalization condition individually satisfied by 
proton and neutron as 

2 / dV / (1^ / Sh{uj,p>,Pp{r)) dio = Z, 

2y / Sh{u:,Tp,pn{r)) dio = A-Z , (A18) 

1 /3 

For a symmetric nuclear matter of density p(r), there is a unique Fermi momentum given by pf{^) = [37r^p(r)/2] 
for which we obtain 

4 J J Shiuj,p,pF{r))du} = p{r) (A19) 

leading to the normalization condition given by 

'ijd^r J J Shiui,p,p{r)) duj = A, (A20) 

where p{r) is the baryon density for the nucleus which is normalized to A and is taken from the lepton nucleus 

scattering experiments. Also we calculate the average kinetic and total nucleon energy given by 

<T >= ^ J d^r J -^^{E{p)-M) J Shipo,P,pir)) dpo , (A21) 

< E >= ^ J d^r J J Sh{po,P,pir)) Podpo, (A22) 

and the binding energy per nucleon given by [i^: 

\EA\=-li<E-M>+^^<T>) (A23) 

Here we have ensured that we retrieve the kinetic energy < T > and the total energy < A > for the nucleon. We 

have tabulated in Table U the kinetic energy per nucleon and the binding energies for nuclei used in the numerical 

calculations. 


[11 E. D. Bloom et al., Phys. Rev. Lett. 23 , 930 (1969); M. Breidenbach et al., Phys. Rev. Lett. 23 , 935 (1969); H. W. Kendall, 
Rev. Mod. Phys. 63 , 597 (1991). 























19 


[2] R. L. Jaffe and M. Soldate, Phys. Lett. B 105, 467 (1981). 

[3] A. Donnachie and P. V. Landshoff, Phys. Lett. B 95, 437 (1980). 

[4] O. Nachtmann, Nucl. Phys. B 187, 141 (1981). 

[5] G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977), V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972) 
[Yad. Fiz. 15, 781 (1972)], L. N. Lipatov, Sov. J. Nucl. Phys. 20, 94 (1975) [Yad. Fiz. 20, 181 (1974)], Y. L. Dokshitzer, 
Sov. Phys. JETP 46, 641 (1977) (Zh. Eksp. Teor. Fiz. 73, 1216 (1977)]. 

(6] J. J. Aubert et a7]The European Muon Gollaboration], Phys. Lett. B 123, 275 1983. 

(7] R. G. Arnold et ai, Phys. Rev. Lett. 52, 727 (1984). 

(8] A. G. Benvenuti et aZ.jThe BGDMS Collaboration], Phys. Lett. B 189, 483 (1987). 

(9] D. M. Aide et ai, Phys. Rev. Lett. 64, 2479 (1990). 

(10] J. Ashman et aZ.]The European Muon Collaboration], Z. Phys. C 57, 211 (1993). 

(11] J. Gomez et ai, Phys. Rev. D 49, 4348 (1994). 

]12] P. Amaudruz et a/.jThe New Muon Collaboration], Nucl. Phys. B 441, 3 (1995). 

(13] K. Ackerstaff et al. (HERMES Collaboration], Phys. Lett. B 475, 386 (2000) (Erratum-ibid. B 567, 339 (2003)] 
arXiv:hep-ex/9910071 . 

(14] L. H. Tao et ai, Z. Phys. C 70, 387 (1996). 

]15] M. Arneodo et aZ.(The New Muon Collaboration], Nucl. Phys. B 487, 3 (1997). 

(16] M. E. Christy et al. (E94110 Collaboration], Phys. Rev. C 70, 015206 (2004). 

(17] V. Tvaskis et al, Phys. Rev. Lett. 98, 142301 (2007). 

]18] V. Mamyan, arXiv: 1202.1457 (nucl-ex]. 

(19] P. Solvignon, D. Gaskell and J. Arrington, AIP Conf. Proc. 1160, 155 (2009) arXiv:0906.0512 ]nucl-ex]]. 

]20] J. Seely, A. Daniel, D. Gaskell, J. Arrington, N. Fomin, P. Solvignon, R. Asaturyan and F. Benmokhtar et al., Phys. Rev. 
Lett. 103, 202301 (2009). 

]21] A. Daniel, Ph.D. thesis. University of Houston, 2007. 

]22] P. Amaudruz et al. (New Muon Collaboration], Phys. Lett. B 294, 120 (1992). 

(23] P. Gueye, M. Bernheim, J. F. Danel, J. E. Ducret, L. Lakehal-Ayat, J. M. Le Goff, A. Magnon and C. Marchand et al., 
Phys. Rev. C 60, 044308 (1999). 

(24] S. Dasu et al., Phys. Rev. Lett. 60, 2591 (1988). 

(25] S. Dasu, P. deBarbaro, A. Bodek, H. Harada, M. W. Krasny, K. Lang, E. M. Riordan and L. Andivahis et al., Phys. Rev. 
D 49, 5641 (1994). 

(26] A. C. Benvenuti et al., Phys. Lett. B 237 (1990) 592. 

(27] M. Arneodo et al. Nucl. Phys. B 483 (1997) 3. 

(28] V. Guzey, L. Zhu, C. E. Keppel, M. E. Christy, D. Gaskell, P. Solvignon and A. Accardi, Phys. Rev. C 86, 045201 (2012). 

(29] 1. Schienbein et al., Phys. Rev. D 77, 054013 (2008), 

Phys. Rev. D 80, 094004 (2009); K. Kovarik et al., 

Phys. Rev. Lett. 106, 122301 (2011). 

(30] J. Arrington, A. Daniel, D. Day, N. Fomin, D. Gaskell and P. Solvignon, Phys. Rev. C 86, 065204 (2012). 

(31] G. A. Miller, S. J. Brodsky and M. Karliner, Phys. Lett. B 481, 245 (2000). 

(32] P. Castorina, Phys. Rev. D 72, 097503 (2005). 

(33] J. J. Zhang, J. H. Gao and X. N. Wang, Phys. Rev. D 91, 014026 (2015). 

(34] D. de Florian and R. Sassot, Phys. Rev. D 69, 074028 (2004). 

]35] M. Hirai, S. Kumano and T. H. Nagai, Phys. Rev. C 76, 065207 (2007). 

(36] K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904, 065 (2009). 

(37] D. de Florian, R. Sassot, P. Zurita and M. Stratmann, Phys. Rev. D 85, 074028 (2012). 

(38] K. Kovarik, 1. Schienbein, T. Stavreva, F. 1. Olness, J. Y. Yu, C. Keppel, J. G. Morfin and J. F. Owens, Few Body Syst. 
52, 271 (2012). 

(39] Akulinichev, S.V. et al. Phys. Lett. B 158 485 (1985) 485, ibid Phys. Rev. Lett. 55 (1985) 2239, ibid Phys. Lett. B 234 
(1990) 170. 

(40] G. V. Dunne and A. W. Thomas, Phys. Rev. D 33, 2061 (1986). 

(41] R. P. Bickerstaff and A. W. Thomas, J. Phys. G 15, 1523 (1989). 

(42] C. Ciofi Degli Atti and S. Liuti, Phys. Lett. B 225, 215 (1989). 

(43] S. A. Kulagin, Nucl. Phys. A 500 (1989) 653. 

(44] M. Arneodo, Phys. Rept. 240, 301 (1994). 

(45] D. F. Geesaman, K. Saito and A. W. Thomas, Ann. Rev. Nucl. Part. Sci. 45 (1995) 337. 

(46] O. Hen, D. W. Higinbotham, G. A. Miller, E. Piasetzky and L. B. Weinstein, Int. J. Mod. Phys. E 22, 1330017 (2013). 

(47] G. Piller and W. Weise, Phys. Rept. 330, 1 (2000). 

(48] E. Marco, E. Oset and P. Fernandez de Cordoba, Nucl. Phys. A 611 (1996) 484. 

(49] O. Benhar, V. R. Pandharipande and 1. Sick, Phys. Lett. B 469, 19 (1999); ibid Phys. Lett. B 410, 79 (1997). 

(50] J. R. Smith and G. A. Miller, Phys. Rev. Lett. 91, 212301 (2003) (Erratum-ibid. 98, 099902 (2007)]. 

(51] S. A. Kulagin and R. Petti, Nucl. Phys. A 765 (2006) 126. 

(52] C. Ciofi degli Atti, L. L. Frankfurt, L. P. Kaptari and M. 1. Strikman, Phys. Rev. C 76 (2007) 055206. 

]53] S. A. Kulagin and R. Petti, Phys. Rev. D 76 (2007) 094023. 

(54] G. A. Miller, Eur. Phys. J. A 31 (2007) 578. 

(55] M. Sajjad Athar, S. K. Singh and M. J. Vicente Vacas, Phys. Lett. B 668, 133 (2008). 


20 


[56] S. A. Kulagin and R. Petti, Phys. Rev. C 82 (2010) 054614. 

[57] M. Sajjad Athar, I. Ruiz Simo and M. J. Vicente Vacas, Nucl. Phys. A 857, 29 (2011). 

[58] H. Haider, I. Ruiz Simo, M. Sajjad Athar and M. J. Vicente Vacas, Phys. Rev. C 84, 054610 (2011). 

[59] L. Frankfurt and M. Strikman, Int. J. Mod. Phys. E 21, 1230002 (2012). 

[60] H. Haider, I. Ruiz Simo and M. Sajjad Athar, Phys. Rev. C 85, 055201 (2012). 

[61] H. Haider, I. Ruiz Simo and M. Sajjad Athar, Phys. Rev. C 87, 035502 (2013). 

[62] H. Haider, M. Sajjad Athar, S. K. Singh and I. R. Simo, Nucl. Phys. A 940, 138 (2015). 

[63] S. Malace, D. Gaskell, D. W. Higinbotham and I. Cloet, Int. J. Mod. Phys. E 23, 1430013 (2014). 

[64] M. Ericson and A. W. Thomas, Phys. Lett. B 128, 112 (1983). 

[65] R. P. Bickerstaff and G. A. Miller, Phys. Lett. B 168, 409 (1986). 

[66] E. L. Berger and E. Goester, Ann. Rev. Nucl. Part. Sci. 37 (1987) 463. 

[67] R. L. Jaffe, Phys. Rev. Lett. 50, 228 (1983). 

[68] H. Mineo, W. Bentz, N. Ishii, A. W. Thomas and K. Yazaki, Nucl. Phys. A 735, 482 (2004). 

[69] 1. C. Gloet, W. Bentz and A. W. Thomas, Phys. Rev. Lett. 95, 052302 (2005). 

[70] O. Nachtmann and H. J. Pirner, Z. Phys. C 21, 277 (1984). 

[71] F. E. Close, R. G. Roberts and G. G. Ross, Phys. Lett. B 129, 346 (1983). 

[72] L.L. Erankfurt and M.I. Strikman, Phys. Rept. 160, 235 (1988). 

[73] N. Armesto, J. Phys. G 32 (2006) R367. 

[74] A. C. Benvenuti et al. [BCDMS Collaboration], Phys. Lett. B 195, 91 (1987). 

[75] JLab 12 GeV proposal PR12-11-113, spokespersons: L. Y. Zhu, M. E. Chirsty, C. E. Keppel, D. Gaskell, and P. Solvignon, 
Detailed Studies of the Nuclear Dependence of R = nn/crT . 

[76] https://www.jlab.org/12GeV/updates.html 

[77] P. Eernandez de Cordoba and E. Oset, Phys. Rev. C 46 (1992) 1697. 

[78] C. Garcia-Recio, J. Nieves and E. Oset, Phys. Rev. C 51 (1995) 237. 

[79] 1. Schienbein et al., J. Phys. G 35 (2008) 053101. 

[80] Pavel M. Nadolsky et al., Phys. Rev. D 78, 013004 (2008); http://hep.pa.msu.edu/cteq/public 

[81] J. A. M. Vermaseren et al., Nucl. Phys. B 724 (2005) 3. 

[82] W. L. van Neerven and A. Vogt, Nucl. Phys. B 568 (2000) 263; ibid 588 (2000) 345. 

[83] M. Gluck, E. Reya and A. Vogt, Z. Phys. C 53 (1992) 651. 

[84] M. Gluck, E. Reya and 1. Schienbein, Eur. Phys. J. C 10 (1999) 313. 

[85] C. G. Callan, Jr. and D. J. Gross, Phys. Rev. Lett. 22 (1969) 156. 

[86] P. Renton, Electroweak Interactions: An Introduction to the Physics of Quarks and Leptons, Cambridge University Press, 
Cambridge (1990). 

[87] S. Moch, J. A. M. Vermaseren and A. Vogt, Phys. Lett. B 606 (2005) 123. 

[88] C. Itzykson and J. B. Zuber, Quantum Pield Theory, McGRAW-HILL Publication, Fourth Edition, Singapore(1988). 

[89] C. Mahaux, P. F. Bortignon, R. A. Broglia and C. H. Dasso, Phys. Rept. 120, 1 (1985). 

[90] L. L. Frankfurt and M. 1. Strikman, Phys. Lett. B 183, 254 (1987). 

[91] H. De Vries, C. W. De Jager and C. De Vries, Atom. Data Nucl. Data Tabl. 36 (1987) 495. 

[92] C. Garcia-Recio, J. Nieves and E. Oset, Nucl. Phys. A 547, 473 (1992). 


